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Due to an excellent mass resolution and a large amount of available data, the CDF experiment, 
located at the Tevatron proton-antiproton accelerator, allows the precise measurement of spectro- 
scopic properties, like mass and decay width, of a variety of states. This was exploited to examine 
the first orbital excitations of the baryon, the resonances At(2595) and At(2625), in the de- 
cay channel A+ 7r+ 7t^, as well as the A^ spin excitations Ec(2455) and Ec(2520) in its decays 
to A+ and A+ n'^ final states in a data sample corresponding to an integrated luminosity of 
5.2fb^'. We present measurements of the mass differences with respect to the Ac and the decay 
widths of these states, using significantly higher statistics than previous experiments. 
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1. Introduction 

Heavy quark baryons provide an interesting laboratory for studying and testing Quantum Chro- 
modynamics (QCD), the theory of strong interactions. The heavy quark states test regions of the 
QCD where perturbation calculations cannot be used, and many different approaches to solve the 
theory were developed. In this write-up we focus on Ac (2595), Ac (2625), rc(2455) and Ec(2520) 
baryons. All four states were observed before and some of their properties measured [|l|]. How- 
ever, different mass measurements for rc(2520) and Ac(2625) are inconsistent and the statistics 
for Ac(2595) and Ac(2625) is rather low. In addition, Blechman and co-workers showed that a 
more sophisticated treatment, which would take into account the proximity of the threshold in the 
Ac(2595) decay, yields a Ac(2595) mass which is 2 — 3MeV/c^ lower than the one observed 
Ec states were observed and studied in AcTT decays, while excited Ac states decay mainly to a 
Ac;r;r final state and decays through intermediated Zc resonances are possible. One peculiarity of 
the experimental studies of these baryons is in their cross talks, which requires special care in the 
treatment of the background due to different kinematic regions allowed for different sources. 

In this analysis we exploit a large sample of A+ — >■ pK^n^ decays collected by the CDF 
detector to perform the measurement of the masses and widths of the discussed charmed baryons. 
We take into account all cross-talks and threshold effects expected in the decays under study |^]. 
Throughout this document the use of a specific particle state implies the use of the charge-conjugate 
state as well. 

2. Analysis Strategy 

The employed dataset corresponds to an integrated luminosity of 5.2 fb^^ and was collected 
using the displaced track trigger which requires two charged particles coming from a single vertex 
with transverse momenta larger than 2GeV/c and impact parameters in the plane transverse to the 
beamline in the region of 0. 1 to 1 mm. 

The selection of the candidates is done in two steps. As all final states feature a Ac daughter, 
as first step we perform a Ac selection. In the second step we perform a dedicated selection of the 
four states under study. In each step we use neural networks to distinguish signal from background. 
All neural networks [Q] are trained using data only by means of the g^^lot technique [^]. As we 
use only data for the neural network trainings, for each case we split the sample to two parts (even 
and odd event numbers) and train two networks. Each of them is applied to the complementary 
subsample in order to maintain a selection which is trained on a sample independent from the one 
to which we apply it. 

In order to determine the mass differences relative to the Ac and the decay widths of the six 
studied states, we perform binned maximum likelihood fits of three separate mass difference (AM) 
distributions. The first two are A+7r+ and A+tt" where the states rc(2455)++'*' and rc(2520)++'° 
are studied. The last one is A+TI+TT^ for Ac(2595)+ and Ac(2625)+. In each of the distributions 
we need to parametrize two signals and several background components. 

Both rc(2455)++'" and rc(2520)++''' are described by a nonrelativistic Breit-Wigner function 
convolved with a resolution function. The resolution function is parametrized with three Gaussians 
centered in zero and other parameters derived from simulated events. The mean width of the 



2 



Charm baryon spectroscopy at CDF 



Felix WICK 



500 



CDF Run II preliminary, L = 5.2 fb 

— Data 

— Fit Function 
■■■E^{2455r-i.A; n* 
■■■Ic{2520r^A; n* 

N(Z,(2455)'*) = 13800 
N(Z,(2520)*') = 8800 



150 200 250 300 



a) 



Mass(Ac n^)-Mass(Ac) [MeV/c^] 



g_ 2000 



-g 1500 



CDF Run II preliminary, L = 5.2 fi 

— Data 

— Fil Function 
■■■1^(2455)°^ a; n 
■■■1^(2520)°^ a; n 

N(I,(2455)'') = 15900 
N(I,(2520)'') = 9000 




b) 



Mass(Aj ii>Mass(A^) [MeV/c^] 



Figure 1: The (a) M(A+;r+) -M(A+) and (b) M(A+;r^) — M(A+) distribution obtained from data (points 
with error bars) together with the fit projection (black Hne). The brown and purple lines correspond to the 
two signal contributions, the green line represents the combinatorial background without real Ac, the blue 
line shows real Ac combined with a random pion and the red dotted line represents a reflection from excited 
Ac decays. The red dashed line corresponds to the sum of all three background contributions. 



resolution function is about 1.6MeV/c2 for 1^(2455)++'*' and about 1.6M&N/c^ for 1^(2520)++'". 
We consider three different types of background, namely random combinations without real A^, 
combinations of real A^ with a random pion and events due to the decay of excited A^ to A+tt+tt^. 
Thereby, random combinations without a real A^ dominate and are described by a second-order 
polynomial with shape and normalization derived in a fit to the AM distribution from A^ mass 
sidebands. The apparent difference between doubly-charged and neutral combinations is due to 
D*(2010)+ mesons with multibody decays, where not all D° decay products are reconstructed. 
In order to describe this reflection, an additional Gaussian function is used. The second background 
source, consisting of real A^ combined with a random pion, is modeled by a third-order polynomial. 
As we do not have an independent proxy for this source, all parameters are left free in the fit. The 
last source, originating from nonresonant Ac(2625) decays, is described by the projections of the 
flat A+TT+TT^ Dalitz plot on the appropriate A+7r+ respective A+tt^ axes. The measured data 
distributions together with the fit projections can be found in figure 

In the fit for Ac(2595) and Ac(2625), an additional complication compared to the case 
arises from the fact that previous measurements of the Ac(2595) properties indicate that it decays 
dominantly to ££.(2455)71, which has its threshold very close to the Ac (2595) mass. Ref. [g] shows 
that taking into account the resulting strong variation of the natural width yields a lower Ac(2595) 
mass than observed by former experiments. With CDF statistics we are much more sensitive to the 
details of the Af (2595) line shape than previous analyses. The Ac(2595) parametrization follows 
Ref. [^. The state is described by a nonrelativistic Breit-Wigner function with a mass-dependent 
decay width which is calculated as integral over the Dalitz plot containing the £^(2455) reso- 
nances. The pion coupling constant h2 is related to the Ac (2595) decay width and represents the 
quantity we measure instead of the natural width. The shape is then numerically convolved with a 
three Gaussian resolution function determined from simulation, which has a mean width of about 
1.8MeV/c^. The signal function for the Ac (2625) is a nonrelativistic Breit-Wigner which is again 
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Figure 2: The M(A+;r+;r^) — M(A+) distribution obtained from data (points with error bars) together with 
the fit projection (black line). The brown and purple lines correspond to the two signal contributions, the 
green line represents the combinatorial background without real A^, the blue line shows real Ac combined 
with two random pions and the red dotted line represents real Ec(2455) combined with a random pion. The 
red dashed line corresponds to the sum of all three background contributions. 
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convolved with a three Gaussian resolution function determined from simulation. The mean width 
of the resolution function is about 2.4MeV/c^. The background consists of three different sources 
which are combinatorial background without real A^, real combined with two random pions 
and real SfJ^^'" combined with a random pion. The combinatorial background without real Ac is 
parametrized by a second order polynomial whose parameters are determined in a fit to the AM 
distribution of candidates from A^. mass sidebands. The second source, consisting of real A^. com- 
bined with two random pions, is parametrized by a second order polynomial with all parameters 
allowed to float in the fit. The model of the final background source, real (2455) combined with 
a random pion, is based on a uniform function defined from the threshold to the end of the fit range, 
where the natural widths as well as the resolution effects of the rc(2455)++'° are taken into account 
for the threshold determination. The size of this contribution is constrained to the rc(2455) yield 
obtained from fits to the projections of the A+tt+tt^ Dahtz plot from the upper A^ (2625) mass 
sideband on the appropriate A+7r+ respective A+tt^ axes. The measured data distribution together 
with the fit projection can be found in figure ^. 

We investigate several systematic effects which can affect our measurements. Generally, they 
can be categorized as imperfect modeling by the simulation, imperfect knowledge of the momen- 
tum scale of the detector, ambiguities in the fit model and uncertainties on the external inputs to 
the fit of the Ac(2595) signal. The most critical point is that we need to understand the intrinsic 
resolution of the detector in order to properly describe the signal shapes, because an uncertainty in 
the resolution has a large impact on the determination of the natural widths. As we estimate it using 
simulated events, it is necessary to substantiate that the resolution obtained from simulation agrees 
with the one in real data. We use D*(2010)+ D^K^ with D° -> K^n+ and ^if{2S) J/^fK^K^ 
with J/y ^ for this purpose and compare the resolution in data and simulated events as 

a function of the transverse momentum of the pions added to or J/y- We also compare the 
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overall resolution scale between data and simulated events and find that all variations are below 
20% which we assign as uncertainty on our knowledge of the resolution function. 

3. Results 

We select about 13800 1^(2455)++, 15900 1^(2455)°, 8800 £^(2520)++, 9000 1^(2520)°, 
3500 Ac(2595)+ and 6200 Ac(2625)+ signal events and obtain the following results for the mass 
differences relative to the A^. mass and the decay widths: 

AM(r,,(2455)++) = 167.44 ± 0.04 (stat.) ±0.12 (syst.)MeV/c^ 
r(r,,(2455)++) = 2.34 ±0. 13 (stat.)± 0.45 (syst.)MeV/c^ 
AM(rc(2455)°) = 167.28 ± 0.03 (stat.) ±0.12 (syst.)MeV/c^, 
r(Zc(2455)°) = 1.65 ±0.11 (stat.) ± 0.49 (syst.)MeV/c^ 
AM(Ec(2520)++) = 230.73 ± 0.56 (stat.) ±0.16 (syst.)MeV/c^ 
r(r,,(2520)++) = 15.03±2.12(stat.)±1.36(syst.)MeV/c^ 
AM(rc(2520)°) = 232.88 ± 0.43 (stat.) ±0.16 (syst.)MeV/c^ 

r(rc(2520)°) = 12.51 ± 1.82 (stat.)± 1.37 (syst.)MeV/c^ 
AM(Ac(2595)+) = 305.79 ±0. 14 (stat.) ±0.20 (syst.)MeV/c^, 

/j^(Ac(2595)+) = 0.36 ± 0.04 (stat.)± 0.07 (syst.), 
AM(A,.(2625)+) = 34 1. 65 ± 0.04 (stat. ) ±0.12 (syst.) MeV/c^. 

For the width of the At(2625) we get a value consistent with zero and therefore calculate an upper 
limit using a Bayesian approach with a uniform prior restricted to positive values. At the 90% 
C.L. we obtain r(Ac(2625)+) < 0.97MeV/c^. For easier comparison to previous results, h\ cor- 
responds to r(Ac(2595)+) = 2.59 ± 0.30 ± 0.47 MeV/c^. 

Except of AM(Ac.(2595)+), all measured quantities are in agreement with previous measure- 
ments. For AM (Af. (2595)+) we observe a value which is by 3.1 MeV/c^ smaller than the existing 
world average. This difference is of the same size as found in Ref. [Q]. The precision for the 
states is comparable to the precision of the world averages. In the case of excited A^ states, our 
results provide a significant improvement in precision against previous measurements. 
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